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A Comparison  of  Theoretical  and  Experimental  Results 

in  Supercooled  Stratus  Dispersal 


1.  INTRODUCTION 


[ 


The  dissipation  of  supercooled  stratus  clouds  by  seeding  with  dry  ice  or  silver 
iodide  has  often  been  a haphazard  procedure.  Vickers  and  Church1  determined  an 
optimum  size  and  seeding  rate  for  dry  ice  by  applying  operations  analysis  techniques 
to  a series  of  experiments  in  thin  stratus  clouds.  However,  indiscriminate  extrapo- 
lation of  their  results  to  thicker  cloud  decks  or  to  silver  iodide  seeding  can  produce 
overseeding  of  the  clouds.  At  its  worst,  this  prevents  clearing,  because  the  ice 
crystals  generated  remain  suspended  in  the  cloud  instead  of  precipitating  out.  Even 
when  this  does  not  occur,  the  efficiency  of  the  silver  iodide  seeding  technique  is 
decreased  to  the  point  where  the  acquired  visibilities  are  not  adequate  for  the  de- 
sired operation. 

2 

Chappell  and  Smith  derived  a one-dimensional  mathematical  model  of  the 
growth  of  ice  crystals  in  supercooled  clouds.  The  ideal  ice  concentration  is  a 
function  of  temperature,  liquid  water  content,  drop  size  distribution,  vertical  mo- 
tion, and  cloud  thickness.  The  amount  of  seeding  material  required  to  obtain  this 
ideal  ice  crystal  concentration  depends  on  the  cloud  thickness  and  the  efficiency  of 
the  seeding  material  which  is,  in  turn,  a function  of  temperature. 

(Received  for  publication  9 August  197  8) 

1.  Vickers,  W.W.,  and  Church,  J.  (1966)  Investigation  of  optical  design  for  super- 

cooled cloud  dispersal  equipment  and  techniques,  J.  Meteorol.  5(No.  1):  1 05- 118. 

2.  Chappell,  C.F.,  and  Smith,  D.  R.  (1976)  On  the  Crystal  Concentrations  Needed 

to  Dissipate  Cold  Stratus  Clouds,  NOAA,  Atmospheric  Physics  and  Chemistry 

Laboratory,  Boulder,  Colorado  (Unpublished  report). 
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It  is  usually  not  possible  to  measure  all  pertinent  cloud  parameters  prior  to 
seeding  and  to  date  there  has  been  no  direct  verification  of  the  Chappell -Smith 
model.  However,  seeding  experiments  using  chlorine-doped  silver  iodide  flares, 
conducted  in  Northern  Michigan  in  February  1977,  included  measurements  of  cloud 
depths  and  temperatures  as  well  as  observations  of  clearing  time.  In  this  report, 
a sensitivity  analysis  is  done  with  the  model  to  determine  the  important  parameters 
and  how  they  affect  the  clearing  efficiency.  Model  calculations  are  then  compared 
with  actual  field  results. 


2.  SENSITIVITY  ANALYSIS  OF  THE  MATHEMATICAL  MODEL 

The  equations  and  assumptions  entering  the  Chappell-Smith,  one-dimensional 
time-dependent,  microphysical  model  are  described  in  Appendix  A.  This  appendix 
is  part  of  an  unpublished  report  submitted  by  NOAA/APCL  to  AFGL  in  response  to 
work  done  for  AFGL  under  Project  Order  No.  Y7T-807.  In  this  model,  ice  nuclei 
of  a given  concentration  are  distributed  uniformly  throughout  a cloud  with  a given 
temperature,  liquid  water  content,  drop  concentration,  and  vertical  velocity.  The 
ice  crystals  then  grow  by  diffusion  and  accretion  until  they  fall  out  of  the  cloud  layer. 
Fall  speeds,  based  on  crystal  mass  and  updraft  velocity,  are  calculated  in  order 
to  determine  fallout  times.  A sensitivity  analysis  was  performed  on  this  model  to 
determine  the  effect  of  varying  individual  parameters.  Some  of  the  results  derived 
from  this  analysis  are  presented  here. 

2.1  Seeding  Rate 

Pioneers  in  cloud  dispersal  sometimes  had  the  naive  belief  that,  if  they  were 
able  to  produce  a clearing  after  a long  interval,  a higher  seeding  rate  would  have 
produced  the  same  clearing  more  quickly.  This  would  be  true  if  the  only  governing 
factor  were  the  rate  at  which  the  cloud  drops  are  converted  to  ice  crystals.  However, 
calculations  show  that  the  visibility  restriction  due  to  ice  crystals  is  greater  than 
that  due  to  the  liquid  water.  Therefore,  the  crystals  must  either  fall  out  or  evapor- 
ate below  the  cloud  base  before  clearing  can  take  place.  The  fallout  time  is 
directly  related  to  how  fast  the  ice  crystals  take  on  water— the  fewer  the  crystals, 
the  less  competition  for  the  available  water  and  the  faster  they  will  grow  and  fall  out. 

This  is  illustrated  in  Figure  1.  The  curves  of  this  figure  were  obtained  by 
varying  the  number  of  ice  nuclei  per  volume  of  air,  keeping  all  other  parameters 
constant.  The  time  it  takes  to  achieve  an  arbitrary  horizontal  visibility  of  1.  5 km 
steadily  decreases  with  increasing  nuclei  concentration.  The  resulting  horizontal 
visibility  will  vary  with  depth  because  of  the  differences  in  the  ice  crystal  residence 
time.  If  the  seeding  material  is  distributed  uniformly  over  the  depth  of  the  cloud. 
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then  residence  time  decreases  with  increased  height  and,  therefore,  maximum 
visibility  improvement  occurs  near  the  base  of  the  cloud. 


t 

i 


TIME  (min) 


Figure  1,  Effect  of  Ice 
Nuclei  Concentration  on 
Clearing  Rate 


The  time  it  takes  for  the  resulting  ice  crystals  to  fall  through  the  cloud  actually 
increases  with  increasing  seeding  rate.  As  can  be  seen  in  Figure  1,  increasing  the 
seeding  rate  will  only  increase  this  time.  For  the  thinner  cloud  (250  m)  the  seeding 
rate  is  not  so  critical.  At  seeding  rates  less  than  3 nuclei/ liter,  the  ice  nuclei  fall 
out  of  the  thinner  cloud  before  the  visibility  near  the  cloud  base  reaches  1.  5 km,  thus 
restricting  the  visibility  improvement  to  something  less  than  1.  5 km.  However,  in 
such  shallow  cloud  layers,  the  ground  can  be  visible  even  when  the  cloud  visibility 
is  much  less  than  1.5  km. 


2.2  Seeding  Rate 


Cloud  temperature  affects  the  efficiency  of  the  seeding  material  (the  number  of 
ice  nuclei  generated  per  gram  of  seeding  material)  and  also  the  rate  of  diffusion 
of  water  vapor  to  the  ice  crystals.  The  model  assumes  a given  concentration  of  ice 
nuclei  per  volume  and,  therefore,  treats  only  the  second  effect.  The  diffusion  rate 
is  a function  of  the  difference  between  the  vapor  pressure  over  water  and  over  ice. 
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Since  the  greatest  difference  occurs  at  approximately  -12°C,  it  would  be  expected 
that  the  faster  clearings  would  occur  at  this  temperature.  Figure  2 shows  this  to 
be  true.  However,  the  differences  in  clearing  time  are  quite  small  between  -10° 
and  -25°C.  At  temperatures  warmer  than  -4°C  the  vapor  pressure  difference  is 
so  small  that,  at  least  for  the  given  conditions,  droplets  are  unable  to  evaporate 
sufficiently  to  raise  the  visibility  to  a level  greater  than  1.  5 km.  Again  it  is  noted 
that  the  fallout  time,  and  not  the  evaporation  time,  defines  the  clearing  time. 


TIME  (min) 


Figure  2.  Effect  of  Temperature  on  Clearing  Rate 

2.3  Vertical  Motions 

Computations  at  a variety  of  temperatures,  seeding  rates,  and  liquid  water 
contents  were  made  assuming  updrafts  of  up  to  10  cm/sec.  As  expected,  the  fallout 
time  increases  with  increasing  updraft  velocity,  thus  increasing  the  time  to  clear 
as  shown  in  Figure  3.  In  this  example,  the  clearing  time  increases  23  percent  (from 
44  min  to  54  min)  when  the  updraft  velocity  increases  from  0 to  10  cm/sec. 


I 
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LWC-  2fl/ir.» 
CRYSTAL  CONC  = SO/liler 


Figure  3.  Effect  of  Vertical  Motion 
on  Clearing  Rate 


2.4  Liquid  Water  Content  and  Cloud  Drop  Size 

When  computations  were  made  in  which  the  cloud  liquid  water  content  (LWC)  was 
varied,  the  results  were  as  shown  in  the  solid  curve  of  Figure  4— a gradually  de- 
creasing time  to  clear  with  increasing  LWC.  In  these  computations  the  cloud  drop 
size  was  held  constant,  so  that  increasing  LWC  also  means  Increasing  concentration. 
The  decreasing  time  with  increasing  liquid  water  is  due  to  the  fact  that  the  crystals 
grow  larger  with  increasing  liquid  water  and  thus  have  a greater  terminal  velocity. 
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Figure  4.  Effect  of  Cloud 
Liquid  Water  Content  and 
Drop  Size  on  Clearing  Rate 
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Assuming  no  accretion  or  riming,  cloud  drop  size  has  an  insignificant  effect 
on  clearing  rate.  However,  with  increasing  drop  size  and  liquid  water  content,  the 
probability  of  accretion  and  riming  increases.  Larger  cloud  droplets  have  higher 
collision  efficiencies,  and  the  higher  the  liquid  water  content,  or  concentration,  the 
more  likely  that  collisions  and  riming  will  occur.  For  a given  mass,  rimed  crystals 
fall  faster  than  unrimed  crystals  and,  therefore,  will  fall  out  sooner.  Figure  4 
shows  that  riming  occurs  when  cloud  droplets  are  greater  than  4 gm  radius  and  the 
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liquid  water  content  is  greater  than  about  0.  3 g/m  . For  cloud  droplets  below 

2.  5 gm  radius  the  model  assumes  the  collection  efficiency  to  be  zero.  Between 
2.  5 and  4 gm,  the  collection  efficiencies  are  quite  low  but  not  zero,  thus  causing 
some  accretion  but  not  enough  to  cause  riming— at  least  as  defined  by  the  model. 

The  model  arbitrarily  assumes  that  riming  occurs  if  the  accretion  rate  is  greater 
than  the  diffusion  rate.  No  experimental  data  exist  to  demonstrate  how  accurate  this 
assumption  is  and  whether  or  not  the  sharp  discontinuity  shown  in  Figure  4 actually 
occurs.  In  reality,  there  is  probably  a gradual  transition  from  unrimed  to  rimed 
particles. 


3.  COMPARISON  BETWEEN  MODEL  AND  FIELD  TESTS 

The  tests  conducted  in  northern  Michigan  early  in  1977  have  been  described 
elsewhere.  ’  3  4 In  this  section  we  shall  summarize  the  results,  and  compare  them 
with  predictions  based  on  the  mathematical  model.  The  actual  ice  nuclei  concentra- 
tion, which  is  an  input  parameter  of  the  model,  depends  on  the  amount  of  seeding 
material  dispensed  per  unit  volume  and  its  efficiency,  which  is  a function  of  tempera- 
ture. The  solid  line  in  Figure  5 is  the  efficiency  curve  for  the  chlorine- doped 
silver  iodide  flares  used  in  the  1977  Michigan  program.  The  results  are  expressed 
as  nuclei  per  g-am  of  seeding  material.  Seeding  rates  were  expressed  as  grams 
of  Agl  per  nautical  mile  per  1000  ft  cloud  depth.  Assuming  a seeding  width  of 
1000  m and  converting  to  ice  nuclei  per  liter,  the  units  used  in  the  mathematical 
model,  we  obtain 

nuclei/liter  = nuclei/nautical  mile/  1000_  ft  . (1) 

*.  649  X 10  1 


3.  Dyer,  R.  M.  et  al  (1977)  Dispersal  of  Supercooled  Stratus  Clouds  by  Silver 

Iodide  Seeding,  Preprints,  Sixth  Conference  on  Planned  and  Inadvertent 
Weather  Modification,  Amer.  Meteor.  Soc. , p.  184. 

4.  Wisner,  C.,  and  Thompson,  J.R.  (1978)  Supercooled  Stratus  Dispersal  Test 

Program,  Final  Report,  AFGL  Contract  No.  F19623-76-C-0306. 
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Figure  5.  Efficiency  of 
Chlorine-doped  Silver  Iodide 
Flares  (solid  line!  and 
Non-doped  Flares  (dashed  line) 


Table  1 summarizes  the  field  results  and  compares  them  with  the  model  cal- 
culations. The  "Observed"  parameters  are  those  parameters  actually  measured 
during  the  field  tests.  Cloud  depth  was  determined  from  altitude  measurements  by 
the  aircraft  flying  at  cloud  base  and  at  cloud  top.  The  cloud  top  and  cloud  base 
temperatures  were  obtained  from  a temperature  probe  mounted  on  the  aircraft. 

The  seeding  rate  in  terms  of  g/nautical  mile/ 1000  ft  depth  is  obtained  from  the 
number  of  flares  ejected.  In  the  single  line  tests,  as  compared  with  multiple  line 
tests,  a variable  seeding  rate  was  used.  Conversion  of  this  seeding  rate  to  nuclei 
per  liter  was  done  by  using  the  efficiency  curve  of  Figure  5,  and  the  formula  of 
Eq.  (1).  Time  to  clear  in  the  "Observed"  section  is  the  time  at  which  the  observer 
in  the  observation  plane  reported  seeing  the  ground  through  the  stratus  deck.  The 
main  source  of  error  in  this  section  is  the  conversion  from  seeding  rate  to  crystal 
concentration,  depending  as  it  does  on  a laboratory-measured  efficiency  curve  with 
some  scatter  in  the  measurements.  Another  source  of  imprecision  is  the  observed 
time  to  clear.  The  observation  aircraft  did  not  always  remain  atop  the  seeded  area, 
and  on  those  occasions  the  clearing  would  have  occurred  some  minutes  before  the 
observation. 
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Table  1.  Comparison  Between  Theoretical  and  Observed  Clearing  Times 


Min 

Cloud 

Depth 

Observed 

Seeding  Rate 

Time 

to 

Clear 

Calculated 
Time  to 
Clear 

TempCC) 

(m) 

(g/nmi/kft) 

(No.  ID 

(min) 

(min) 

-4 

400 

20-200 

. 007  -.  07 

No  Hole 

No  Hole 

-6 

270 

20-200 

. 638-6.38 

No  Hole 

No  Hole 

-8.5 

1890 

30-300 

53.  1-531 

Partial^  ** 

45 ->  80 

-8.5 

1070 

30 

53 

39 

27  -51 

-9 

1200 

10-100 

35.  4-354 

Partial 

29-69 

-9.  5 

1370 

80 

354 

70 

42-77 

-9.  5 

1200 

200 

88.  6-886 

>36<2) 

32 ->  80 

-10 

1310 

70 

433 

70 

64  ->  80 

-11.  5 

7 60 

36 

383 

45 

38-44 

-12.  0 

520 

7-70 

99.  2-992 

11(3) 

29-55 

-14.0 

460 

9-90 

287-2870 

30 

16-44 

-15.  5 

730 

8-80 

354-3543 

38 

36 ->  80 

-17.  5 

910 

44 

2728 

38 

>80 

-17.5 

270 

30-300 

1860-18600 

20 

20 ->  80 

-21.  5 

488 

5-50 

708-7086 

45 

27  ->  80 

(1)  20  g flares  were  used.  This  size  flare  burns  over  only  a 1200  m depth. 

(2)  Log  does  not  indicate  when  ground  became  visible,  but  it  was  sometime 
after  36  minutes. 

(3)  Thin  spots  were  observed  in  cloud. 


The  computer  model  requires  values  of  cloud  liquid  water  content,  drop  size, 
and  concentration  and  vertical  motion,  but  these  parameters  were  not  measured. 
Calculations  were  made  using  the  observed  conditions  and  a variety  of  assumed 

liquid  water  content  and  drop  size  conditions.  Liquid  water  content  varied  from 

3 3 

.05  - 0.  5 g/m  and  drop  diameter  from  10.  8-24  jum.  The  0.  5 g/m  liquid  water 

3 

contents  represent  cases  in  which  riming  occurred,  and  the  .05  g/m  represents 
cases  in  which  some  accretion  took  place  but  no  riming,  according  to  the  model. 

3 

A cloud  drop  concentration  of  7 5 /cm  and  a zero  vertical  velocity  were  assumed. 

As  shown  earlier,  a change  in  the  drop  concentration  has  little  impact  on  the  clear- 
ing time  except  where  it  means  the  difference  between  riming  and  no  riming.  Also 
cases  with  updrafts  will  increase  the  clearing  time  slightly. 

The  calculated  clearing  times  show  a wide  range  of  possible  times.  Computa- 
tions were  stopped  at  80  min,  so  the  notation  >80  indicates  that  clearing  had  not 


occurred  up  to  the  time  computation  stopped.  In  general,  the  slowest  clearing 
occurs  at  low  liquid  water  content  and  high  seeding  rate.  In  eight  of  the  eleven 
cases  in  which  the  ground  became  visible,  the  observed  clearing  times  were  within 
the  calculated  range.  The  model  successfully  predicted  the  two  cases  in  which  no 
clearing  was  observed.  Although  the  seeding  rate  is  relatively  high,  the  number  of 
nuclei  produced  at  these  temperatures,  as  illustrated  in  Figure  2,  is  quite  low. 

The  model  shows  the  ice  crystals  falling  out  before  the  droplets  can  evaporate  enough 
to  raise  the  visibility  sufficiently  to  see  through  the  cloud  layer. 

At  the  lower  temperatures  (<  -15°C)  and  at  the  higher  seeding  rates  the  model 
predicts  overseeding  with  ice  crystals  remaining  in  suspension  rather  than  precipi- 
tating out.  This  may  account  for  the  annulus  of  unprecipitated  ice  crystals  noted  in 
the  center  of  the  seeded  area  where  crystals  concentration  is  highest.  It  is  more 
likely  that  the  actual  ice  crystal  concentration  near  the  edges  of  the  affected  area 
were  considerably  less  than  near  the  center,  thus  resulting  in  more  rapid  clearing 
and  higher  visibilities  along  the  edge.  The  model,  being  a one -dimensional  model, 
does  not  take  into  account  horizontal  diffusion  and  assumes  a uniform  distribution 
of  ice  crystals. 

4.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  field  tests  have  demonstrated  the  feasibility  of  silver  iodide  seeding  for 
tactical  dispersal  of  supercooled  stratus  clouds,  and  the  one -dimensional  mathe- 
matical model  has  provided  reasonable  estimates  of  the  clearing  rate  under  a 
variety  of  meteorological  conditions. 

Further  field  tests,  including  measurements  of  cloud  microphysics,  are  re- 
quired to  derive  a feasible  method  of  determining  optimum  seeding  procedures 
under  operational  conditions.  A field  program  using  a tactical  aircraft  would  be 
the  ultimate  test  of  whether  this  technique  can  become  a standard  Air  Force 
capability. 

The  mathematical  model  needs  extension  and  refinement.  Results  of  field  tests 
during  which  cloud  physics  parameters  are  measured  can  be  used  to  refine  the 
equations  used  in  the  model.  Extension  of  the  model  to  two -dimensions  will  allow 
prediction  of  the  ultimate  size  of  the  opening,  and  perhaps  how  long  the  visibility 
remains  above  minimum.  Incorporation  of  viewing  angle  into  the  model  would  also 
be  desirable,  so  that  slant  range  visibility  may  be  computed. 

Laboratory  studies  and  field  tests  should  also  be  directed  toward  determining 
the  optimum  composition  of  the  silver  iodide  flares.  The  flares  used  in  the  Air  Force 
tests  in  Michigan  contained  a chlorine  compound,  which  improved  their  efficiency. 

The  difference  between  the  solid  and  the  dashed  curve  of  Figure  5 is  the  difference 
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between  the  efficiencies  of  the  chlorine -treated  and  the  untreated  flares.  There 
are  many  unanswered  questions,  such  as  what  proportion  of  chlorine  will  optimize 
the  efficiency  and  what  mechanism  is  operating  to  produce  the  ice  nuclei  from  the 
chlorine-doped  silver  iodide. 
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Appendix  A 


Mathematical  Modal  of  Crystal  Growth  in  Clouds 


Al.  INTRODUCTION 

2 

The  model  derived  by  Chappell  and  Smith  has  not  been  described  in  the  open 
literature.  A brief  summary  of  the  equations  used,  taken  from  Chappell  and  Smith' 
internal  report  to  AFGL,  is  presented  here. 


A 2.  LIST  OF  SYMBOLS 

a Droplet  curvature  term  = 3.  3(10  ^)/T 

b Droplet  solubility  term  = 1.  33925(10 

C Specific  heat  of  air  at  constant  pressure 

Cv  Specific  heat  of  air  at  constant  volume 

D Diffusivity  of  water  vapor 

Dj  Diameter  of  water  droplet 

D Diameter  of  partially  rimed  crystal 

P 

D Diameter  of  unrimed  crystal 

u 

e Vapor  pressure  of  the  environment 

E Collection  efficiency  of  crystals  for  supercooled  cloud  water 

e . e Saturation  vapor  pressure  of  the  environment  (over  a plane 

l w 

ice  surface,  over  a plane  water  surface) 
a Gravitational  acceleration 
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n 


G 

Gi 

K 


Lc*Ls 

mc'md 

N ,N. 
c d 


cw 


R ..R 
a v 

S..S 
l w 

t 

T 

V 

V_ 


W 

t 


P 

♦ i*  wv% 

♦l**2‘*3 


Thermodynamic  function  in  droplet  growth  equation 
Thermodynamic  function  in  ice  crystal  growth  equation 
Thermal  conductivity  of  air  or  scattering  efficiency  of  a 
cloud  droplet 

Latent  heat  of  condensation  and  sublimation 

Mass  of  ice  crystal  and  water  droplet 

Concentration  of  ice  crystals  and  water  droplets 

per  unit  volume  of  air 

Atmospheric  pressure 

Cloud  liquid  water  per  unit  volume 

Radius  of  water  droplet 

Specific  gas  constant  for  dry  air,  and  for  water  vapor 
Supersaturation  with  respect  to  ice  and  water 
Time 

Temperature 

Visibility 

Fall  velocity  of  partially  rimed  ice  crystals 
Fall  velocity  of  unrimed  ice  crystals 
Updraft  velocity 

Ratio  of  the  molecular  weight  of  water  vapor  to  that 
of  air  (0.622) 

Density  of  air 

Thermodynamic  functions  in  the  supersaturation  equation 
Thermodynamic  functions  in  the  supersaturation  equation 


1 

U 


! 


I 


A3.  CLOUD  SUPERSATURATION 

Temporal  changes  in  the  supersaturation  of  a mixed-phase  cloud  are  given  by 

Sj  = OjW  - 02(NdI"d  + VV  • ♦3(Nc”c  + Kme} 

- (S.  + 1)  (1  - P4Ndmd  - <?5Ncmc)  (p/p)  <A1) 

where  m_,  and  m are  the  time  rate  of  changes  in  droplet  and  crystal  mass,  and 
• .d  c 

and  N are  the  time  rate  of  changes  in  the  droplet  and  crystal  concentrations, 
d c 

N is  assumed  equal  to  zero,  since  agglomeration  is  not  considered.  Dots  indicate 
differentiation  with  respect  to  time.  Eq.  (Al)  is  obtained  by  combining  the  Clausius- 
Clapeyron  relation,  energy  equation,  and  the  equations  of  state  for  dry  and  moist 
air.  It  describes  changes  in  cloud  supersaturation  as  a function  of  the  interaction 
among  the  updraft,  crystal  growth,  and  droplet  growth  or  evaporation. 
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The  thermodynamic  functions  in  Eq.  (Al)  are  given  by 


»2-(e^r  *if-)  • 

. "dT\ 

'■'w  r)' 

/ 6 L L \ 


■w- 


The  last  parenthetical  expression  in  Eq.  (Al)  is  given  by 

- [*  iw  + +2{Nd™d  + W ^a^c  + Ncmc)] 
P/P-  1 -*2Ndmd  '^3Ncmc  ’ 


where  the  thermodynamic  functions  are  expressed  by 

C gw 
i v6 

*1  " F/T  * 

d p 
L.  R , 

*2=  C^' 


L R. 

\1/  - s d 
C p ‘ 
P 


A4.  DROPLET  GROWTH 


Droplet  growth  in  Eq.  (Al)  is  evaluated  with  the  expression 


md=4nrdG  Sw  - (a/r^  + (b/rd>  . 


_ c 

The  Kelvin  curvature  term  is  evaluated  by  a = 3.  3(10  )/T,  and  the  nucleus 

solubility  term  is  b = 1.  33925(10  The  thermodynamic  function  is  given  by 

G = |(L^/KRvT2)+  (RvT/ewD)|  . (A12) 

AS.  CRYSTAL  GROWTH 

For  unrimed  or  slightly  rimed  crystals,  depositional  and  accretional  growth 
equations  are 

m , = 205.  195m1/2S.G.(l  + 24.  63m0,  305)  (A13) 

cd  c l l c ' ' 


and 


m = 1709.  26 Eq  V m , 
ca  c w u c 


(A14) 


where  ice  crystals  are  assumed  to  be  planar  disks  and  ventilation  effects  are  in- 
cluded. The  thermodynamic  function  is 


Gi  = (L2/KRvT2)  + (RvT/eiD) 


(A15) 


and  the  crystal  fall  velocity  is 


V 


u 


62U0'  217 
u 


(A16) 


The  relationship  between  crystal  size  and  mass  incorporated  into  Eq.  (A  13)  is 
mc  = 0.  00038D2  . (A17) 


The  collection  efficiency  of  crystals  for  supercooled  cloud  droplets  is  defined 
by  the  function 

E = 0 for  Dd<5p  (A18) 

E = -0,  58576916  + 0.  13872401Dd  - 0.  0030517645D2  for  5p  * Dd  < 17.  5p 
E = 0.  91  for  Dd>  17.  5p  . 
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Crystal  growth  in  Eq.  (A  1>  must  include  only  ice  growth  by  deposition.  How- 
ever, the  crystal  mass  in  Eq.  (A13)  is  modified  by  the  accretional  growth  rate  given 
by  Eq.  (A  14).  The  time  rate  of  change  of  droplet  concentration  due  to  accretion 


is  then 


N . = - N m /m  . . (A19) 

d c ca  d 

For  partially  rimed  crystals  the  depositional  and  accretional  growth  equations 


m . = 96.  225m1(,2S.G.<l  + 28.  57  m0-  325) 
cd  c i i c 


m = 240.  56  Eq  V m . 
ca  Mcw  p c 


The  crystal  fall  velocity  is  given  by 


0 ^ 

V = 17  0.  7 Du‘  J 
P P 


and  the  relationship  between  crystal  size  and  mass  incorporated  into  Eq.  (A20)  is 


m = 0.  0027D  . 
c p 


A6.  VISIBILITY 


Visibility  due  to  the  presence  of  cloud  droplets  is  computed  from 
N 

V = 3.  912/Tr  D (Km.r2  ) , 
i=i  i 


where  N is  the  number  concentration  of  cloud  droplets,  r^  is  the  cloud  droplet 
radius,  and  K is  the  scattering  efficiency  of  a cloud  droplet.  The  scattering  efficiency 
is  taken  to  have  a value  of  approximately  two  for  visible  light  and  spherical  cloud 
droplets. 


